Critical voltage of a mesoscopic superconductor 
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We study the role of the quasiparticle distribution function / on the properties of a superconduct- 
ing nanowire. We employ a numerical calculation based upon the Usadel equation. Going beyond 
linear response, we find a non-thermal distribution for / caused by an applied bias voltage. We 
demonstrate that the even part of / (the energy mode /l) drives a first order transition from the 
superconducting state to the normal state irrespective of the current. 
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The energy distribution function of quasiparticles in a 
normal metal is under equilibrium conditions given by 
the Fermi-Dirac distribution /q. In recent years it has 
been demonstrated that in a voltage (F)-biased meso- 
scopic wire (length L) a two-step non-equilibrium distri- 
bution develops 1] with additional rounding by quasi- 
particle scattering due to spin- flip and/or Goulomb in- 
teractions 2]. Figure^ shows the distribution, which 
resembles two shifted Fermi-Dirac functions: 

/(x, e) = (l- x)fo{e + eV/2) + xfo{e - eV/2) (1) 

with e the quasiparticle energy and x the coordinate 
along the wire. For strong enough relaxation {L L^f)^ 
with L(f) is the phase coherence length) and/or high tem- 
peratures {ksT ^ eV) the distribution returns to a 
Fermi-Dirac distribution with a local effective temper- 
ature. 

If the normal wire is replaced by a superconducting 
wire, the attractive interaction between electrons leads 
to the superconducting state. The questions we address 
here are how the distribution function is modified (for a 
typical result see Fig. QJd) and how this affects observable 
properties such as the current-voltage characteristics of 
the system and the breakdown of the superconducting 
state. To relate the distribution function to observable 
quantities, it is convenient to separate the symmetric part 
/l (energy mode) from the asymmetric part fr (charge 
mode) which each have a different spatial and spectral 
form (Fig. and d). In particular we will show that 
the breakdown is characterized by a voltage rather than 
by a current; in other words, the system cannot be triv- 
ially treated as two resistors modelling the normal- to 
supercurrent conversion, with a superconducting element 
characterized by its depairing current in-between. 

The transport and spectral properties of dirty super- 
conducting systems {£e <C ^o, with £e the elastic mean 
free path and the superconducting phase coherence 
length) are described by the quasiclassical Green func- 
tions obeying the Usadel equation Q. For out of equi- 
librium systems we use the Keldysh technique in Nambu 
(particle-hole) space. We look at s-wave superconduc- 



tors (singlet pairing) without any spin-dependent in- 
teractions. The Usadel equation then takes the form 
hDV{GVG) = -i[H,G], where the check notation (G) 
denotes a 4 x 4 matrix, D is the diffusion constant, V is 
the spatial derivative 4] and we neglect any inelastic pro- 
cess. The elements of G and when split up in Keldysh 
space, are 2x2 matrices in Nambu space, denoted by a 
hat: 
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Here, G^ and G^ are the retarded and advanced com- 
ponents describing equilibrium properties and G^ is the 
Keldysh component which describes the non-equilibrium 
properties. Their elements are the quasiclassical (energy- 
dependent) normal and anomalous Green functions and, 
for the Keldysh component only, the quasiparticle dis- 
tribution functions (which take account of the non- 
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FIG. 1: Quasiparticle distribution function f{x,e) as func- 
tion of energy e and position x for a normal wire (a) and a 
superconducting wire (b) between normal metallic reservoirs 
for ksT <^ eV < Ao, with (c) and (d) the decomposition of 
(b) into the charge mode fr and energy mode /l. 
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equilibrium). For the Hamiltonian H we write: 

H={ - 1 (3) 



e 

A* 



where e is the (eigen)energy and the chosen gauge is such 
that the pair potential A is in equilibrium a real quan- 
tity, A = A*. The matrix Green function G satisfies the 
normalization condition GG = 1, leading to G^G^ = 
G^G^ = 1 and G^G^ + G^G^ = 6. If supercon- 
ducting reservoirs in the system are kept at zero voltage 
(avoiding AC Josephson effects), G^ can be written as 
G^ = G^ f — fG^ . Here / is the diagonal generalized dis- 
tribution number matrix of the quasiparticles in Nambu 
space. To relate / to observable quantities we decompose 
it into an even part (or energy/longitudinal mode) and 
an odd part (or charge/transverse mode) in particle- hole 
space: / = /l^o + /t^s, where are the Pauli matrices 
in particle- hole space 5]. The full distribution function 
is retained by: 2f{x^£) = 1 — /l(^,^) — frix^s). 

The retarded matrix Green function in terms of the po- 
sition and energy dependent normal g{s^x) and anoma- 
lous Fi{e^x) Green functions is: 



qR ^ ( di^^x) Fi{e,x) 
F2{e,x) -g{e,x) 



(4) 



Substituting this in the retarded part of the Usadel equa- 
tion: hDV(G^VG^) = -i[H, G^] and using the normal- 
ization condition (^^ + F1F2 = 1), we find the retarded 
Usadel equations: 



hD[g\/^Fi - FiV^g] = -2iAg - 2ieFi 
hD[FiV'^F2- F2\/'^Fi] = 2iAF2 + 2iA*Fi 



(5) 



The second equation is essential when calculating the 
non-equilibrium properties of superconductors. Its left- 
hand-side is proportional to the divergence of the spectral 
(energy-dependent) supercurrent, which is (compared to 
the equilibrium case) no longer a conserved quantity. 
A general relation between the advanced matrix Green 
function and the retarded matrix Green function is given 
by: = —T3{G^)^T3. Using this, the Keldysh matrix 
Green function can be written entirely in terms of g, 
Fi, F2, fL and /r: 
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(6) 



where f± = ^ Jt- Working out the kinetic part 
of the Usadel equation: hDV{G^VG^ + G^VG^) = 
—i[H^ G^] we find (combining the diagonal components) 
the kinetic equations describing the non-equilibrium 
part: 



hDVjenergy = 

hDVjch^rge = "^RlJl + 2RTfT 



The various elements in Eq. [3 are given by: 



Jenergy 
Jcharge 

Hi 
Ht 

je 

Rl 

Rt 



niV/L + nxV/T+Je/T 

n^v/T-nxV/i+je/i 

2\g\^-\F,\'-m^) 



■2i5r + ii^ii 



\F2?) 



= \{\F,\'-\F,\') 
= isR{FiVF2 -F2VF1} 



(8) 



i5{AF2 + AF/} 



Equations [7| are two coupled diffusion equations for 
and /t, describing the divergences in the spectral energy 
current and the spectral charge current. The total charge 
current is given by J = J jcharge<^^ with p the resistiv- 
ity. The terms Ul and IIt can be related to an effective 
diffusion constant for the energy and charge mode respec- 
tively and Ux as a "cross-diffusion" between them, is 
the spectral supercurrent and Rl and Rt describe the 
"leakage" of spectral current to different energies, where 
the total leakage-current (X J[RLfL-\-RTfT]de is zero. In 
the small signal limit the terms Ux , je and Rl are small 
and can in many cases be neglected (linear approach), 
effectively decoupling /l and fr- In this article we go 
beyond this limit. 

The Usadel equation is supplemented by a self- 
consistency relation: 



-H(l,2) - — / <^(l,2)«e 



(9) 



Here, A^o is the normal density of states around the Fermi 
energy, V^f^ the effective attractive interaction and the in- 
tegral limits are set by the Debye energy hwD- The result- 
ing equation for A becomes: A = — ^A^oV^eff j-^^ [(^1 ~ 
Fl)fL-{F,+Fl)fT]de. 

To calculate spectral and transport properties, one 
needs to know the self-consistent solution of A. In most 
practical cases, this has to be done numerically. A con- 
venient solution scheme is to first find the Green func- 
tions of the system by solving the retarded equations 
for a certain A, next to determine the quasiparticle dis- 
tribution functions by solving the kinetic equations and 
then calculate a new A using the self-consistency rela- 
tion. This process has to be repeated until A converges. 
As a starting value for A we use the BCS form at zero 
temperature. To simplify the calculations a parameter- 
ization is used that automatically fulfills the normaliza- 
tion condition. It is convenient to take g = cosh(^), 
Fi = sinh(6>)e*^ and F2 = — sinh(6>)e~*^, where and 
X are position and energy dependent (complex) vari- 
ables. At the interfaces between the superconducting 
wire and the normal metallic reservoirs we use the fol- 
lowing boundary conditions: = Vx = (retarded 
equation) and /l,t = |(tanh ±tanh |^^) (kinetic 
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FIG. 2: The calculated current (J)-voltage(V) relation of a 
superconducting wire of length L = 8.5^o between normal 
metallic reservoirs (see inset) at several temperatures, and 
for a wire of length 17^o at T = 0. Jc is the critical current 
density, and Ao the bulk gap energy 



equation) , where the latter are the usual reservoir distri- 
bution functions. 

The transport properties of the NSN system can now 
be calculated with the equations described above. In 
a previous analysis a finite differential conductance was 
found at zero bias employing a linear response calculation 
0. With the approach introduced here, the full current- 
voltage relation can be obtained. The result at several 
temperatures is displayed in Fig. [SJ with the voltage 
normalized to Ao(= Abuik,T=o) and the current density 
normalized to the critical current density Jc ^ 0.75^^ 

, with ^0 = At T = we observe a linear re- 

sistance at low voltages caused by the decay of /t (Fig. 
d), and a critical point (voltage) above which the resis- 
tance is equal to the normal state resistance. At higher 
temperatures (T = 0.5,0.75Tc) a linear approach would 
only give an adequate approximation in a limited volt- 
age range. We will argue below that the superconductor 
switches to the normal state by which is controlled 
by the voltage and cannot be interpreted as a critical 
current. 

In Fig. El the electrostatic potential (j) = fT^{g}de 
along the wire is shown at zero temperature prior 
to (eV/Ao = 0.013,0.646) and immediately after 
(eV/Ao = 0.651) the transition. The potential can be 
seen to drop to zero over a distance of the order of the 
coherence length due to the normal- to supercurrent con- 
version. This mechanism also gives rise to the finite zero 
bias resistance. The profile hardly changes over the full 
range of voltages, until the critical value is reached, after 
which the electrostatic potential drops in a linear fash- 
ion, indicating the system is in the normal state. The 
minimal changes emphasize the limited influence of fr 



on the superconducting state (i.e. on A). 

The current density at which the superconductor 
switches to the normal state (for T = 0) is much smaller 
than the critical current density in an infinitely long wire 
(J/Jc = 1). Neither is the transition triggered at the 
weaker superconducting edges as indicated by the shape 
of the electrostatic potential profile in Fig. |S1 

The parameter that determines whether or not the su- 
perconducting state exist is A, as follows from Eq. El 
The integral in this self-consistency equation sums all 
pair-states (either occupied by a Cooper pair, or empty). 
Fi gives the Cooper pair density-of-states and /l and 
/t determine which of those states are doubly occupied 
or doubly empty and which are singly occupied (broken) 
due to the presence of quasiparticles. In equilibrium at 
T = 0, a switch to the normal state can only be caused 
by reaching a critical phase gradient, entering A via F^. 
In the presence of quasiparticles, A (and thus potentially 
the state of the system), is also influenced by the distri- 
bution functions. It was noticed above that the charge 
mode /t has a very limited influence on A. The effect of 
the energy mode /l is examined below. 

By a small modification of our system to a T-shaped 
geometry as shown in Fig. ^ we can in a direct way dis- 
entangle the effects of /l and fr on A. This setup can 
be thought of as the connection of the superconducting 
wire to the center of a normal wire. In the middle of 
such a wire fr is equal to zero, but /l is not. The result 
for the pair potential at the edge of the superconduct- 
ing wire as a function of the voltage of the reservoirs is 
shown in Fig. ^ Although there is no net current flow- 
ing through the superconductor, at a certain voltage the 
pair potential collapses. The voltage that is necessary to 
trigger this transition to the normal state is very close to 
the transition in Fig. [21 (where we used the two termi- 
nal setup). Apparently the influence of /l is important, 
since it can cause the superconductor to switch to the 
normal state irrespective of the value of the supercurrent. 
Clearly the influence of /l on the state of the supercon- 
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FIG. 3: The normalized electrostatic potential as a function 
of position x along the superconducting wire for bias voltages 
prior to and immediately after the transition (at T = 0) . 
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FIG. 4: Top: T-shaped geometry, bottom: pair potential 
A in the S-wire as function of (left) voltage at two differ- 
ent positions; (right) position at two different voltages. The 
breakdown voltage is at eV/ Aq = 0.707 

ductor is larger than the influence of the supercurrent on 
this same quantity. 

The quantity that defines the possible states of the 
system is the free energy. Evidently the superconductor 
compares two states for the minimization of this free en- 
ergy: the first state is the superconducting state in which 
the free energy remains constant as a function of voltage 
(and independent of the shape of /l provided this shape 
does not change for energies larger than the gap). The 
second possible state is the normal state. At zero tem- 
perature, in the absence of a bias voltage, the difference 
in free energy between the two states is the condensa- 
tion energy of the superconductor. When the voltage is 
increased (but still eV < A), the free energy of the super- 
conducting state remains constant while the free energy 
of the normal state decreases since in that case electrons 
occupy higher energy states due to the applied voltage. 

To illustrate the effect, we calculate explicitly the free 
energy difference between the superconducting state (Fs) 
and normal state (Fn) at zero temperature for a bulk su- 
perconductor (analytically) and the T-shaped structure 
(numerically), as a function of voltage (which appears in 
/l) and A. At zero temperature the free energy of the 
system reduces to the internal energy (kinetic plus po- 
tential) 01 • From the analytical calculation for the bulk 
following Bardeen |^ we find that /l changes the free 
energy in such a way that at eV = |a/2Ao the supercon- 
ductor undergoes a first order phase transition. For the 
voltage range ^Aq < eV < Aq the state of the system 
has two solutions. The free energy difference for the bulk 
superconductor is shown in Fig. [51|9|. Numerical results 
for the free energy of the T-shaped geometry are shown 
as well (both as function of position and as function of 
voltage). For long wires, the numerical results approach 
the analytical (bulk) calculation. This indicates that the 
effect of the bias voltage can indeed be related to the 
existence of a first order phase transition at zero temper- 
ature. 



T-shape (numerical) bulk (analytical) 




FIG. 5: Free energy difference between the superconducting 
and normal state. Right: analytical bulk solution showing 
the bistable voltage range. Left: numerical solutions for (top) 
increasing wire length as function of voltage and (bottom) as 
function of position. Free energies are normalized to (0)/87r 

In conclusion, we have studied the role of the energy 
mode /l of the quasiparticle distribution on the prop- 
erties of a superconducting nanowire. We employ a nu- 
merical simulation of the Usadel equation in full-response 
and find a non-thermal distribution for /l (caused by an 
applied bias voltage) which drives a first order transi- 
tion from the superconducting state to the normal state 
irrespective of the current. A direct calculation on the 
free energy of a bulk superconductor confirms that the 
voltage indeed causes the phase transition. In general, 
the significant role played by /l found in these super- 
conducting nanowires stresses the importance of treating 
/l and /t on equal footing. 
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